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ABSTRACT 

The ROSAT soft X-ray spectrum and pulse profile of the Geminga pulsar are analyzed 
and interpreted in terms of thermal emission from the surface of the neutron star. The X- 
ray spectrum appears to consist of two blackbody components with T\ = (5.2 ± 1.0) x 10 s 
K and T2 ~ 3 X 10 6 K, respectively. The inferred ratio of surface areas, A2/A1, is ~ 3 X 
10 -6 . Both components are highly modulated at the pulsar rotation period, but the harder 
X-ray pulse is narrower, and leads the main (soft) X-ray pulse by about 105° of phase. The 
soft X-ray component is interpreted as photospheric cooling of much of the neutron star’s 
surface area, while the small, hot region could be part of the much smaller polar cap heated 
by energetic particles flowing inward from the magnetospheric accelerator which is responsible 
for the production of Geminga’s 7-rays. Geminga’s 7-ray emission is consistent with outer- 
magnetosphere accelerator models for highly inclined dipoles. These predict the beaming of 
energetic 7-rays close enough to the star to give copious e* production in the stellar magnetic 
field and a large circumstellar pair density from pair inflow toward the surface. These pairs 
may quench radio emission, and also reflect most of the hard polar cap X-rays back to the 
stellar surface by cyclotron resonance scattering. They are then reemitted from that much 
larger area at the lower temperature T\. The single-peaked nature of the X-ray pulse and its 
energy-dependent phase suggest an off-center dipole geometry for the surface magnetic field. 

Under the assumption that the soft X-ray emission comes from the full surface of a neutron 
star of radius R = 10 km, a distance estimate of (150-400) pc is derived. This range is consistent 
with the fit interstellar column density of (1.5± 0.5) x 10 20 cm -2 . Distances less than 150 pc are 
probably ruled out both by the lower limit on the column density, and also by the requirement 
that the Rayleigh- Jeans extrapolation of the soft X-ray spectrum not exceed the observed blue 
flux of the faint optical counterpart. This distance estimate implies that Geminga’s efficiency 
for converting spindown power into 7-rays is near unity, and that there may be significant 
beaming of the 7-rays as well. These results tend to bolster the prospect that most of the 
unidentified high-energy 7-ray sources in the Galactic plane are pulsars, some of which may be 
radio quiet. 
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1. INTRODUCTION 

The high-energy 7-ray source Geminga (2CG 195+04, IE 0630+178) is a radio-quiet pulsar 
with period P = 0.237 s, surface field B p ~ 1.6 X 10 12 G, and characteristic age r = P/2P = 
3.4 x 10 5 yr (Halpern & Holt 1992; Bertsch et al. 1992). Its soft X-ray properties contain a great 
deal of information about thermal emission from the surface of the neutron star. There are very 
few pulsars for which the bulk of the soft X-ray emission is demonstrably photospheric, and it 
is argued here that Geminga is one of the more secure cases. Therefor* , it can be used to test 
theories of neutron star cooling, polar-cap heating, and the effect of the surface magnetic field 
and environment on the emergent spectrum. Other pulsars for which i uch surface emission is 
likely to contribute substantially are PSR 0656+14 (Finley, Ogelman, < 1 Kiziloglu 1992), PSR 
1055-52 (Cheng & Helfand 1983; Brinkmann & Ogelman 1987), and V«jla (Ogelman, Finley, & 
Zimmermann 1992), but Geminga allows the most detailed examination so far of the spectrum 
and pulse properties of the photospheric emission. 

If interpreted as blackbody emission, the X-ray spectrum also provides the only measure- 
ment of the distance to Geminga, which is crucial to the determination of the efficiency for 
turning spindown power into 7-rays. This in turn affects the likelihood that the other uniden- 
tified Galactic high-energy 7-ray sources are pulsars. Geminga bolsters the hypothesis that the 
COS B sources are a population of efficient 7-ray pulsars. The X-ray pulse profile of Geminga 
is complex and highly modulated, and to the extent that it can be interpreted in terms of the 
geometry of the magnetic field, it may help in understanding the 7-ray emission mechanism 
and pulse profile. This paper presents a more comprehensive analysis of the ROSAT observa- 
tion of the Geminga pulsar than was possible in the initial report by Halpern fc Holt (1992). 
Although the data are the same, a number of additional results have been obtained, and the 
interpretation also incorporates properties of the 7-ray emission. 

2. PULSE PROFILE AND TIMING 

A more detailed soft X-ray pulse profile was constructed by folding about the revised 
EGRET ephemeris (P and P) of Mayer-Hasselwander et al. (1992), which is more precise 
than the original X-ray determination. Table 1 gives the folding parameters which apply to the 
epoch of the ROSAT observation in the Barycentric Dynamical Time (TDB) system. A total of 
7,911 counts in the 0.07-1.50 keV band, of which less than 2% are background, were extracted 
from a circle of radius 150" in the position sensitive proportional counter. Figure 1 shows a 
number of details that were not apparent in the original examination of the X-ray data. First, 
there is apparently a sharp dip in the soft X-ray flux, and possibly significant high-frequency 
structure, especially on the rising side of the pulse. A revised estimate of the pulsed fraction 
in the 0.07-0.28 keV band is 33%. Although it appeared in Halpern and Holt (1992) that the 
pulsed fraction decreases with energy, it is now seen that there is actuaUy a faint but highly 
modulated pulse at energies around 1 keV which leads the peak of the main X-ray pulse by 
about 105° of phase. The weakness of the modulation around 0.5 keV is due to the fortuitous 
superposition of these two out-of-phase components, and it appears that both components are 
individually highly pulsed in a roughly energy-independent manner. We shaU refer to the ~ 0.2 
keV X-rays and the ~ 1 keV X-rays as the “soft” and “hard” components respectively. The 
pulsed fraction of the hard X-rays is ~ 37%, similar to that of the soft X-ray pulse, but the 
hard pulse is narrower. The dip at phase 0.5 is not apparent in the hard light curve, although 
the photon statistics here are quite poor. The spectral analysis presented in the next section 
lends additional support to the interpretation of the pulse profiles as two separate components. 
In particular, the soft component will be interpreted as emission from the full surface of the 
neutron star, while the hard X-ray pulse could arise from the heated polar caps. 

The relative phasing of the X-ray and 7-ray pulses is of great interest, since the pulse 
profiles are so different. The best current analysis of the ROSAT spacecraft clock leads to 
the alignment shown in Figure 2. The time of phase 0.0 is given by To in Table 1, which is 
an extrapolation of the EGRET ephemeris over the five- week interval between the ROSAT 



observation and the first EGRET observation. The phase comparison using P and P should be 
accurate to better than 0.01 cycles in the absence of significant timing noise, and is comparable 
to the accuracy of the ROSAT spacecraft clock (Ogelman, Finley, & Zimmermann 1992). [It 
should be noted, however, that P as given by Hermsen et al. (1992) and Mattox et al. (1992) 
probably is dominated by timing noise or cycle count errors, since the quoted values lead to 
unreasonably large values of the braking index, 31 and 44, respectively.] Another possible 
source of uncertainty is the position of Geminga, which could cause an error of < 7 msec (0.03 
cycles) in the arrival time at the barycenter if it is off by < 3* in right ascension. An error 
of this magnitude would be expected if the G" optical counterpart (Halpern & Tytler 1988) is 
incorrect, or if there is exceptionally high proper motion. The tentative detection of a proper 
motion of ~ 0."2 yr -1 by Bignami, Caraveo, & Mereghetti (1992) would seem to alleviate both 
of these concerns. More worrisome is a possible ambiguity of 0.5 s in the interpretation of the 
times from the ROSAT spacecraft clock (Ogelman 1992, private communication), which would 
result in an error of 0.11 cycles in the phase of the pulse. For now, we tentatively accept the 
relative phasing of the X-ray and 7-ray pulses shown in Figure 2 as correct to within ~ 0.13 
cycles. 

The sharp, double-peaked 7-ray pulse is in marked contrast to the broad, single-peaked 
X-ray pulse. In § 5.2, it is argued that Geminga’s high 7-ray efficiency and the 180° separation 
of its 7-ray peaks call for a highly inclined dipole. If so, the single hard X-ray pulse cannot 
be understood in terms of polar-cap emission from a dipole close to the center of the star, 
since both polar caps should be visible with a near 180° separation of their X-ray emission. 
The soft X-ray component is presumed to come from nearly the full surface of the star, and it 
too is single-peaked. The simplest resolution of this problem invokes an off-center dipole, as 
described in § 5.5, so that the polar caps are too close together to be resolved in the light curve 
of the rotating star. The 7-ray pulses originating in the outer magnetosphere can nevertheless 
be double fan beams which are visible from most directions. The hard X-ray pulse roughly 
coincides with the bridge of 7-ray emission between the two main peaks, which is consistent 
with an off-center dipole configuration if the 7-rays are coming from regions along the open 
field lines which are roughly symmetric with respect to the polar caps. 

3. SPECTRAL ANALYSIS 

The X-ray source counts were extracted from a circle of radius 150", and background from 
an annulus between radii of 150* and 250*. A total of 7,637 net counts were detected in the 
0.08-2.45 keV band in a live time of 14,205 s. As shown by Halpern & Holt (1992), the spectrum 
cannot be fitted by a single blackbody or power-law, both of which leave a residual hard excess 
above 0.6 keV. Two-component models are adequate, and one which is particularly amenable 
to interpretation is a pair of blackbodies, 

F(E) = C E 3 + eE/ ^ — ) keVcm" 2 s" 1 keV 1 . (1) 

The interstellar medium photoelectric cross sections of Morrison & McCammon (1983) axe 
used for <r(E). The four fitting parameters are the temperatures T\ and Tj* interstellar column 
density Ah, and fraction / which is closely related to the ratio of surface areas (A 2 /Ai) of the 
two blackbody components. (For each trial spectrum, the normalization constant C is fixed by 
the total counts, and its error is largely systematic.) A full four-dimensional x* grid search was 
performed, and the best fit is' shown in Figure 3. 

Before interpreting these results, it is necessary to point out that there are a number bf 
sharp features in the residuals between 0.2 and 0.4 keV which have been seen in other sources, 
and which are thought to reflect an as yet unsolved problem in the pulse-height analysis or the 
detector response matrix. These features are not all likely to be real, and they contribute to an 
unacceptable total x 2 of 61 f° r 42 degrees of freedom. In the hope that the overall continuum 



shape is not significantly affected by this problem, a simple renormalization of the x 2 values by 
the factor 42/61 was made to permit sensible confidence contours to be derived. 

A number of interesting results can be drawn from this fit. First, the relative contributions 
of the two blackbodies cross over at ~ 0.6 keV, which is very close to the energy of minimum 
modulation as seen in the light curve (Figure 1). This coincidence lends additional credence 
to the idea that two components are present in the X-rays, since they can be separated both 
spectrally and temporally. The lower temperature blackbody with 7\ ~ 5 X 10 s K accounts for 
about 96% of the bolometric flux, and its fit is relatively insensitive to the presence of the weaker 
hard component. If the x 2 grid is projected onto the (Ti,Nh) plane, the confidence contours 
shown in Figure 4 restrict the temperature and column density to the ranges (3 - 6) X 10 5 K 
and (1 — 3) x 10 20 cm -2 , respectively. The confidence contours are elongated because the errors 
on T\ and Ah are correlated. Further restrictions on these parameters are made possible by 
the incorporation of optical data as described in § 4 below. 

A second interesting projection of the x 2 grid is onto the ( X-j , / ) plane. Figure 5 shows 
that T 2 lies in the range (2.3 - 3.9) X 10 6 K, and that the fractional area / is ~ 3 X 10 -5 . 
The very large range allowed for /,10 -6 - 10~ 4 , is dominated not by the uncertainty in the 
flux of the hard component, but rather by the ranges in Ah, Tj, and normalization of the soft 
component since a large and uncertain fraction of the soft emission is absorbed in the intervening 
interstellar medium. Projections of the corresponding best fit values of 7j are shown as dashed 
lines in Figure 5. A summary of the best fit spectral parameters is given in Table 2. The hard 
component is interpreted in § 5.3 in terms of polar-cap heating of an off-center dipole by an 
outer magnetosphere accelerator. A division of the pulse height spectrum into several time- 
resolved segments according to rotation phase demonstrates that each of the two components 
varies in an approximately energy independent manner. That is, no variations in 7j or T 2 
are detectable as a function of rotation phase, although the paucity of counts above 0.5 keV 
precludes a meaningful test for variation in T 2 . 

A second not implausible two-component model of the X-ray spectrum consists of a black- 
body plus a power law, 

F(E) = C + / S"*) e-°W N » keVcm- 2 s" 1 keV\ (2) 

where a is the energy index. The best fit spectrum is shown in Figure 6, and the parameters 
and uncertainties are listed in Table 3. The fit is nearly as good (x 2 = 63) as the double 
blackbody. Since the power-law component makes a small contribution to the total X-ray 
flux, the fit parameters for the soft blackbody are not significantly different from those in the 
double blackbody model. Although the best fit 7\ is changed slightly, the confidence contours 
encompass the same range of parameters, as shown in Figure 7. An extrapolation of the power- 
law component with F{E) oc A -1-5 to the visible would exceed the flux of the 25 mag G" optical 
counterpart (Halpern & Tytler 1988) by a factor of several hundred. If the hard component 
were to be interpreted as synchrotron emission from the magnetosphere, the local magnetic 
field must be such that the low-frequency cutoff at u = eBfmc lies above the optical band. 
This implies that B > 10 9 G. Therefore, any X-ray synchrotron emission would have to arise 
from within 10 stellar radii of the surface of the neutron star. Synchrotron self-absorption is 
not an alternative explanation for the weakness of the optical emission, since i/*b* < 10* Hz for 
this component. Although formally acceptable, a synchrotron origin for the hard X-rays seems 
less attractive than surface emission because it is difficult to understand why magnetospheric 
emission far from the stellar surface should have a single-peaked pulse. Regardless of the nature 
of the hard X-rays, the weak optical flux could be synchrotron emission coming from radii larger 
than 10 7 cm. The optical flux in the visible band is only 1.7 x 10 -15 ergs cm -2 s -1 , which is a 
factor of 100 less than that of the “hard” X-ray component at 1 keV. 


4. DISTANCE TO GEMINGA 



The softer X-ray component currently provides the only measurement from which the 
distance to Geminga can be inferred. Several constraints on the distance are associated with 
the spectral fit. For each point in Figure 4, the normalization constant C can be converted 
to a visual magnitude by extrapolation of the Rayleigh-Jeans tail of the X-ray spectral fit. If 
we assume in addition that the X-ray emission comes from the full surface of a neutron star 
of radius R = 10 km, then a distance can be assigned to each fit. The straight lines in Figure 
4 indicate these predicted optical magnitudes and distances. Although all the information in 
Figure 4 is derived purely from the X-ray spectrum, a number of additional constraints can 
be brought to bear to further restrict the allowed range of parameters. First, the G" optical 
counterpart has magnitudes V — 25.2 i 0.3 (Halpem & Tytler 1988) and B = 26.5 i 0.5 
(Bignami, Caraveo, & Paul 1988). Since the B - V color is redder than Rayleigh-Jeans, it 
would be unreasonable for the extrapolated X-ray spectrum to exceed the optical brightness. 
This would seem to rule out distances smaller than 150 pc in Figure 4. 

Additional restrictions on the distance can be gleaned from the measured interstellar col- 
umn densities in the direction of Geminga. These come principally from absorption-line studies 
of nearby stars; the derived column densities as a function of distance in the Galactic plane 
have been summarized by Paresce (1984) and Frisch & York (1983). The column density to- 
ward Geminga is small, but it is poorly determined because very few stars have been studied, 
and this direction grazes the edge of the local bubble which is devoid of neutral material to 
large distances in the third Galactic quadrant. What can be most reliably extracted from these 
surveys (in particular, from Figure 1 of Frisch & York and Figure 4 of Paresce) is that Ah does 
not rise to 1 x 10 2 ° cm -2 until a distance of at least 150 pc is reached. It is also possible that a 
distance of ~ 500 pc is necessary for Ah to reach this value. So the allowed range of distances 
can conservatively be restricted to 150-450 pc by comparing the X-ray spectral fit in Figure 4 
with these published measurements of the local interstellar column density. 

A final restriction on the distance comes from the requirement that the apparent 7-ray 
luminosity not exceed the spindown power /OO, or at least not greatly, since some beaming is 
possible. The implied upper limit is ~ 400 pc (Bertsch et al. 1992). So when all the constraints 
are combined, the best estimate of the distance to Geminga is 250 pc, with an allowed range 
of 150-400 pc. The corresponding range of X-ray spectral parameters in the restricted space of 
Figure 4 is T x = (5.2 ± 1.0) X 10 5 K, and A H = (1.5 ± 0.5) x 10 20 cm" 2 . 

Of the assumptions that could affect this distance estimate, one might at first suspect 
the “full surface” condition, since if the assumed area of the emitting surface is incorrect, the 
distances indicated in Figure 4 are changed. But in fact, the acceptable range of distances is 
insensitive to changes in the effective emitting area. Although the individual distance contours 
indicated in Figure 4 would change, the extrapolated optical magnitudes cannot. Neither would 
the fit Ah, the observed interstellar column densities, nor the 7-ray efficiency limit. Similarly, 
the assumed form of the “hard” X-ray component does not significantly affect the fit to the soft 
component, so the distance limits indicated in Figure 7 are the same as those in Figure 4. The 
combination of conditions restricting the distance to between 150 and 400 pc is quite robust, 
and a much more accurate determination of the X-ray spectral parameters would have to be 
made before the uncertainty in distance could be reduced. A better determination of the actual 
run of column density versus distance in the direction toward Geminga would also be helpful 
in pinning down the distance. If the distance estimate presented here is flawed, it is more likely 
to be a result of gross errors in the detector response matrix, or substantial deviations of the 
spectrum from a blackbody as considered in §5.4, neither of which are obviously present. 

5. INTERPRETATIONS 
5.1. Into the Valley of Death 

Si mila.r it.ipii among the 7-ray spectra and most of the light curves of the known 7-ray pulsars 
(Crab, Vela, PSR 1706-44, Geminga) suggest that a similar, powerful accelerator is operating 
in each of their magnetospheres. If this is indeed the case, such an accelerator would have to 
be located in their outer magnetospheres (r >> R ~ 10 6 cm). This would almost certainly be 



true for the source of the Crab pulsar’s optical radiation. Synchrotron emission seems to be the 
only (incoherent) mechanism which is powerful enough to produce the Crab’s optical light, and 
the radiating e* pairs must move relativistically to give the sharply peaked light curve. This 
limits the local magnetic field to B < 10 8 G, which is found only near the light cylinder, at 
r k = c/0 10 8 cm. But the arrival phase of each of the Crab’s two optical pulses is identical 

to those of its X-rays and 7-rays. This strongly implies a common geometrical location for all 
of the radiation very far from the stellar surface. 

Independent arguments for an outer-magnetosphere accelerator come from the COS B 
and EGRET observations of Geminga. First, 7-ray emission remains strong up to 5 GeV 
(Grenier, Hermsen, & Hote 1991). Such 7-rays would be converted to e* pairs if they crossed 
a B > 2 x 10 8 G, and thus would not be expected to escape if produced within much less 
than 10 R of the stellar surface. In addition, the maximum potential drop along B which can 
be sustained within that distance to the star is about 10 12 V. At higher voltages, curvature 
radiated 7-rays would produce so many e^ pairs as they crossed the local B > 10 9 G, that such 
pair production would quench the accelerating potential drop. The maximum possible particle 
flow (No) through any accelerator along the open field-line bundle between Geminga’s surface 
polar cap and its light cylinder is limited to 

No < fi2 ^ pJ * 3 ~ 8 X 10 31 s“\ (3) 

ec ; 

where B p is the dipole field component at the stellar surface. (A current larger than eN 0 
would itself cause a local B near the light cylinder which is larger than that from the star.) 
Therefore, the maximum power from an inner magnetosphere accelerator is eNo X 10 12 V ~ 10 32 
ergs s -1 . But Geminga’s 7-ray flux in the 50-5000 MeV range is ~ 3.1 X 10“ 9 ergs cm “ 2 s -1 
(Grenier et al. 1991), corresponding to a luminosity of 2.3 X 10 34 ergs s -1 if radiation is emitted 
isotropically at the estimated distance of 250 pc. Indeed, any magnetosphere accelerator which 
manages to mobilize a large fraction of a pulsar’s total spindown power ( Jftft ~ 3.3 x 10 34 ergs 
s -1 for Geminga) must be very far from that pulsar’s surface unless the pulsar is near its radio 
death line (Ruderman & Sutherland 1975; Chen & Ruderman 1993). The maximum possible 
potential drop between the surface and the light cylinder of Geminga is 
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2 x 10 14 V, 


( 4 ) 


and almost all of it must be used to account for Geminga’s very large apparent L. T . Because 
the apparent is so close to the theoretical maximum (~ /ftft), it is likely that there is a 
favorable beaming geometry which enhances the flux in oiir direction. 

The “outer-gap” accelerator model was proposed to explain needed pair production and 
associated energetic photon emission in the outer magnetospheres of the Crab and Vela pulsars. 
Such models may have application to all of the 7-ray pulsars. A turn-off for strong 7-ray emission 
from a Vela-like pulsar was predicted at a period P ~ 0.13 s by Ruderman & Cheng (1988). 
This period depended upon the characteristic magnetic field through the accelerator which 
extends along the closed field-line boundary from the “null surface” where ft ■ B = 0 out to the 
light cylinder. For an aligned pulsar, the outer end of the accelerator is at r max = rj c = c/ft. 
The inner end is at r m j„ = (2/3) r) c , but it moves in sharply with increasing angle between the 
stellar dipole moment p and the spin ft until r m j„ ~ (4/9) (p- ft) 2 | px ft | -2 rj c . The turn-off 
period P ~ 0.13 s comes from taking the characteristic B in the accelerator as that at ri c . Chen 
and Ruderman (1993) used B evaluated at r = 0.5 ri c which gives a 7-ray death line defined by 
5 log B p — 12 log P = 69.5. Geminga, with B p = 1.6 x 10 12 G, falls just below this line, i.e., just 

beyond turn-off. However, because a substantial inclination angle between p and ft moves the 
inner end of the accelerator nearer to the 6tar, it increases the accelerator’s magnetic field. This 



could account for Geminga being a strong 7-ray pulsar rather than a moribund one. Indeed, 
if Vela’s inclination angle, 6 = cos -1 (/i • fi), is ~ 35°, and if Geminga’s is ~ 65°, then their 
“outer-gap” accelerators could have similar local magnetic field, size, current flow, and total 
power (Ruderman et al. 1993). Geminga’s accelerator would, however, have to span almost 
all of that star’s accessible open field lines to accomplish this, compared to less than 1/3 of 
the open field lines for Vela’s. With this modified inclination angle, an outer-gap model which 
worked for Vela would probably also be applicable to Geminga. Reports of the death of such 
models (Bertsch et al. 1992) seem exaggerated, or at least premature. 

5.2. Beam Geometry 


Certain features of a Geminga-like outer magnetosphere accelerator seem robust enough to 
survive in any model of an accelerator and 7-ray emission in that region. To motivate these we 
shall, however, use the description of an outer-gap, Vela-like accelerator (Cheng et al. 1986a). 

a) The accelerator, to avoid self-quenching by its own e* pair production, is bounded by the 
last closed field-line surface of Figure 8. There is a similar accelerator on either side of the 
star. 

b) Electrons and positrons accelerated through the accelerator are radiation-reaction limited, 
so that each accelerator is almost 100% efficient as a 7-ray emitter. 

c) To achieve the high L^/ISISI of Geminga, almost all of the appropriately curved open field 
lines which pass through the ft • B = 0 surface also pass through the outer-magnetosphere 
accelerator. 

d) Because of e* pair production within the accelerators, each produces a particle flux N of 
e - (e + ) moving down toward the stair’s polar cap from its inner end, and an equal flux of 
e+ (e~ ) moving outward from its light cylinder. In addition, there is a comparable flux of 
positively (negatively) charged particles pulled up from the polar caps which also passes 
through the accelerator. The total (maximum) flux of e“(e + ) directed down toward the 
polar caps is 

"~ 2 (^) (f£) Ki ~ 4x1031 (5) 


e) 


f) 


g) 


These e“(e + ) will emit the curvature radiation beams 2 and 2' of Figure 8. 

Within the accelerator, pair-produced e - and e + flow in both directions so that the similar 
fan beams of radiation, beams 1 and V from one accelerator, will be matched by beams 
3 and 3' from the other. Because of Geminga’s large tilt angle 6, the curvat are radiation 
beam 2 will only be seen by observers looking within an angle 0 of the pulsai spin axis. It 
is less clear over how large an angle the fan beams 1 and 3 would be observa >le. 

An observer could, in principle, see the three beams 1, 2, and 3 (or F, 2', and 3 1 ). Because of 
Geminga’s large tilt angle, the emission radius r« << ric and time-of-flight a{d aberration 
differences among the three beams axe small. The phase separation between beams 1 and 2 
is A<f> = 0.5 -M with 6 ~ 2r e /xri c . Beams 1 and 3 will merge with separation S. Geminga’s 
observed A<f> and the widths of its two pulses require r t < 0.05 r| C = 5 x 10 7 cm. This 
would imply an inclination angle for Geminga greater than 70°. 

Because of this large angle, beams 2 and 2', and even 1' and 3, may come close enough to 
the strong inner magnetic field of the star to give a very substantial e* pair production 
there. 


5.3. Polar-Cap Heating 

The extreme relativistic particles from the starward end of the accelerator will curvature 
radiate away much of their energy before reaching the stellar polar caps. Their instantaneous 
rate of such energy loss is 


E = 71TIC 2 



( 6 ) 


OF POOH QUALITY 


where r c ~ (rc/ft) 1/2 is the local radius of curvature of the (dipole) field lines which link the 
accelerator to the star’s polar cap. For an accelerator which begins at the radius r min > > R) 
the residual energy of the e^ impacting the polar cap can be approximated as 


E(R) = 



~ 6.5 ergs. 


(7) 


as long as £(r m j„) >> E(R). Then from equation (5), and choosing r min ~ r« ~ 5 x 10 cm, 
the power brought down to the polar cap by particles coming from the inner ends of the outer- 
magnetosphere accelerator is L p ~ 2.6 x 10 32 ergs s -1 . This predicted polar-cap luminosity is 
insensitive to the geometry, location, or pair-production mechanisms in the accelerator, but it 
does assume that the accelerator current (and power) are near maximal. 

If the two heated half polar caps could be directly observed, they would have a temperature 
of ~ 6 X 10 6 K. This is about twice the observed temperature T 2 = (3 ± 1) X 10 6 K of the hotter 
blackbody component. Moreover, the theoretical polar-cap heating exceeds the hard X-ray 
luminosity by a factor of 10 2 , and is even a few times larger than the estimated total X-ray 
emission (see Table 4). While relatively small changes in Geminga’s parameters (e.g., increasing 
the distance by a factor of 2) could accommodate an X-ray luminosity of ~ 2 X 10 , how the 
polar cap radiation is degraded in temperature by a factor of 10 remains a problem. Moreover, 
the observed X-rays seem to be radiated at 5x 10* K from almost the entire neutron star surface 
area of 10 13 cm 2 , instead of from the 10 9 cm 2 polar cap. 

One possible cause for this discrepancy arises in the large density of e* pairs near the 
neutron star. Their effectiveness in reflecting keV X-rays from a hot polar cap is enormously in- 
creased whenever their cyclotron resonance frequency in the strong local magnetic field matches 
the frequency of incident X-rays. We consider next some aspects of this problem. 

5.4. Circumstellar X-ray Scattering 

In the geometry of Figures 8 and 9, the radiation beams 2 and 2' can pass relatively close 
to the neutron star before leaving the magnetosphere. (This may also be true for beams 3 and 
1'.) A 10 GeV 7-ray would be converted into an e* pair by Geminga’s magnetic field if it 
came within 3 X 10 7 cm of the star; a 100 MeV 7-ray would convert if it approached within 
6 x 10 6 cm ~ 6 R. The e ± pairs materialized from the 7-rays of energy E y would themselves 

radiate 7-rays with energy ~ 10(JB^/GeV) MeV. This second generation of 7- rays could 
in turn als o make pairs if some passed within a radius r ~ R(E^ /MeV) 1 / 3 of the star. The 
pair flux thus generated increases with larger inclination angles between Ji and ft, and might 
exceed 10 38 s~* for Geminga. These e* pairs are created mainly on closed field lines and will, 
therefore, ultimately be channeled into the near magnetosphere of the neutron star. Their 
abundance there depends upon the detailed location of the annihilation process, e.g., to what 
extent polar-cap X-ray emission and cyclotron resonant e ± scattering of those X-rays keeps 
these pairs from flowing easily down onto the stellar surface. Local number densities fijt ^ 10 
cm -3 and column densities in excess of 10 21 cm -2 seem quite plausible for Geminga. 

With only a Thomson scattering cross section <r T , these e± would have a negligible optical 
depth to escaping polar-cap X-rays. However, because of the stellar magnetic field the effective 
X-ray e^ scattering cross section can be represented as 


a = B ) 2 + 


27r 2 e 2 

me 


€ x B | 2 6{ub ~ w). 


( 8 ) 


Here E is the (electric field) polarization of the X-ray, and u> B = eB(r)/mc. With fun ~ 0.5 
keV for Geminga’s polar-cap X-rays, and € • B = 0, 




2 x 10 -13 cm 3 . 


(9) 


Thus an n± > 10 13 cm -3 at r ~ 3R would make an optically thick cyclotron resonant backscat- 
tering layer there. Most of the hot polar cap X-rays would then njot escape before being in- 
tercepted by the star. In this way, polar cap X-ray power is transferred to the entire stellar 
surface from which it is reradiated at a predicted lower temperature T \ , where 


r, ~ f di urT * = (£) T; - °i 2Ti 


The predicted ratio Ti/T 2 is in reasonable agreement with the observed values and their errors 
as listed in Table 2. Of course, the cooler surface radiation will have to pass through its own 
cyclotron resonant e^ backscattering layer at r ~ 6 12 before ultimately escaping. However, 
there is no significant further reduction in emitted X-ray energy even if some of this softer 
X-ray flux is also ultimately reabsorbed by the stellar surface, although the spectral shape of 
the emission can be affected. 

A small fraction of the hard polar cap X-rays may pass through the resonant scattering 
layer at r ~ 312 without being backscattered. These are the photons with | t x B |~ 0. They 
are moving almost perpendicular to the local B with polarization along B, so that any cyclotron 
resonances are not effective in amplifying their cross section. Thus, the hard X-ray component 
which passes through the cyclotron resonance backscattering layer might be almost completely 
linearly polarized. 

5.5. X-ray Light Curves 

Any model which tries to explain as other than a coincidence the 180° separation of 
Geminga’s two 7-ray peaks would need a highly inclined dipole. Then, if the stellar sur- 
face magnetic field configuration were nearly that of the canonical central dipole (two surface 
polar caps 180° apart), the hard ~ 1 keV polar cap radiation would be expected to show two 
similar peaks separated by 0.5 in phase. Rather, the data comprise a single peak, or perhaps 
two p eaks sufficiently close together in phase to be unresolvable in data of the present quality. 
This suggests a sunspot-like configuration, such as that from a dipole much nearer the surface 
than a centrally located one. Radiation from a sunspot polar-cap configuration should appear 
brightest to an observer when the star is rotated by about 90° from the phase angle at which 
7-ray beam 2 is observed (Figure 9a). With the highly inclined dipojle and sunspot-like geome- 
try of Figure 9b, this emission direction is also near that for which some of the beam can have 
| € X B |»v 0 (i.e., B perpendicular to the direction of propagation of the photon), and thus 
pass through the circumstellar e* cyclotron resonant scattering layer instead of being scattered 
back toward the star. The maximum hard X-ray intensity would then be expected at a phase 
interval A <j> ~ 0.25 before or after 7-ray beam 2, depending on the sense of the stellar spin. 

The most difficult place for the stellar X-rays to pass through the circumstellar e ± resonance 
layer at 3J2 < r < 6i2 are the areas near the two open field-line bundles, because there B and 
r are nearly parallel. Any X-ray photon moving nearly radially outward there would have 
1 1 X B | *v 1 so that the dominant resonant part of the scattering cross section in equation (8) 
is a maximum. This would cause a minimum in the observed X-ray intensity approximately 
coincident in phase with the y-rays of beam 2 since these come to us from this same open field- 
line bundle direction. The indicated phases of the total X-ray minimum, hard X-ray maximum, 
and 7-ray beam 2 are sketched in Figure 10. 

There is not yet a compelling determination of which of the two observed 7-ray beams 
should be identified with beam 2. However, beam 2 of Figure 8 should consist mainly of 
curvature radiation from e~(e + ) after they have left the accelerator. It would be expected to 
be less intense than the sum of beam 1, which contains similar radiation from e + (e") moving 
outward within and above the accelerator, plus beam 3 from inward moving e (e"*") within 
and above the accelerator on the opposite side of the neutron star. This suggests identifying 
the weaker of Geminga’s two 7-ray pulses (in Figure 10) as beam 2. More compelling support 
would be a softer spectrum at the highest energies for the weaker pulse, as its most energetic 


7-rays are converted to e ± pairs when passing through the stellar magnetic field. If the second, 
stronger 7-ray pulse is mainly beam 3 of Figure 8, then its preceeding interpulse shoulder is 
also suggestive of what is expected from the effects of aberration on emission at larger radii 
(Cheng et ah 1986b). The model phases for Geminga’s X-ray maxima and minima shown in 
Figure 10 are qualitatively similar to the observed ones in Figure 2. 

[A comparison of the expected polar-cap heating with the observed X-ray luminosity can 
also be made for the Vela pulsar. Here the f l and B are larger, but the fraction of open field 
lines subtended by the accelerator is probably of order 10% so that the expected polar-cap 
heating is ~ 4 x 10 32 ergs s“ x , about the same as for Geminga. The observed X-ray luminosity 
in the Vela point source is 4.3 X 10 32 ergs s -x (Ogelman et al. 1992), which is in excellent 
agreement with such a model, but the inferred emitting area of 0.1 times AtB? is more difficult 
to understand: though it is much larger than the polar cap, it is still much smaller than the 
full surface.] 

5.6. Neutron Star Cooling Curves 

In the above interpretation, the soft X-ray component of luminosity ~ 5.4 x 10 31 ergs 
s“ x comes from the stellar surface. Since the observed temperature of 5 X 10 5 K falls on the 
theoretical neutron star cooling curves, there may be no necessary reason to attribute the soft X- 
ray luminosity exclusively (or even mainly) to polar-cap heating. In view of the highly efficient 
conversion of the 3 X 10 34 ergs s _1 spindown power to radiation, it is remarkable that the soft 
X-rays from Geminga still fall within the range allowed by the cooling curves. The predictions 
of a variety of models as presented by Page & Applegate (1992) fall in the range (1 — 7) X 10 5 K 
for a neutron star of age 3 x 10 5 yr, and the Geminga data do not strongly discriminate among 
these. However, Page (1993) concludes that Geminga must have a superfluid core in order for 
its surface temperature to be compatible with most cooling scenarios. Given the expectations 
for polar-cap heating, and the fact that the soft X-ray luminosity is only ~ 10~ 3 of the spindown 
power, the observed temperature T x = (5.2 ± 1.0) X 10 5 K must be regarded as only an upper 
limit when comparing with theoretical cooling curves. An additional reason to interpret T\ as 
an upper limit is the effect of the surface magnetic field on the emergent spectrum, as discussed 
in the next section. 

The shape of Geminga’s soft X-ray pulse is very similar to that of PSR 0656+14 (Finley 
et al. 1992) although the statistics of the latter are poorer because the pulsed fraction is only 
14%. The blackbody temperature of PSR 0656+14 is 9 X 10 s K, and its characteristic age is 
1.1 X 10 5 yr, a factor of 3 younger than Geminga. There is so far no evidence for 7-ray emission 
from PSR 0656+14 ( d 500 pc), nor for an additional hard X-ray component. These facts 

together are consistent with the interpretation of the soft X-ray spectrum in both of these stars 
as “normal” photospheric emission. 

5.7. Magnetized Neutron Star Atmospheres 

The X-ray spectrum emergent from a realistic neutron star atmosphere can deviate sub- 
stantially from a blackbody because of the non-grey opacity (Romani 1987). The observed 
spectrum is generally harder than a blackbody at the given effective temperature. The surface 
magnetic field ameliorates this effect to some extent, but also introduces rotational modulation 
because the radiative opacity and heat conduction are anisotropic (Miller 1992; Shibanov et al. 
1992a, b; Ventura et al. 1992). These authors have modelled both fully ionized and partially 
ionized atmospheres, and the deviations are more substantial for the latter with light element 
composition, since the K-edges of hydrogen and helium are shifted into the ROSAT bandpass 
for typical surface magnetic fields. Although a detailed fitting of these models is beyond the 
scope of this paper, a first-order concern is whether or not magnetic effects can account for the 
complex spectral and pulse behavior in terms of a single component spectrum, without the need 
for two distinct blackbodies as invoked in § 3. The light-element models are the mo6t relevant 
in this regard. 

Beginning with a fully ionized hydrogen photosphere, Shibanov et al. (1992a) showed 
that two models with effective temperature 5 x 10 5 K and B = 4.7 x 10 12 G can differ by as 



much as 45% in flux at 1 keV, according to whether the magnetic field is normal or tangential 
to the surface* At lower energies the difference is less, and below 0.25 keV the spectra cross 
over. Although this anisotropy could in principle be responsible for some of the rotational 
modulation, it seems unlikely to be able to explain the existence of the hard pulse as a tail 
of the soft blackbody, since the observed flux at 1 keV is at least 100 times larger than the 
extrapolation of the soft blackbody component (Figure 3). 

The inclusion of bound-free opacity causes sharper deviations from a blackbody (Miller 
1992; Shibanov et al. 1992b). The theoretical spectrum closely follows a blackbody below the 
photoionization threshold of hydrogen, which is at ~ 200 eV for B ~ 2 x 10 12 G. Deviations can 
be substantial above this energy. In particular, the spectra fall below the blackbody between 
0.2 and 0.6 keV, but rise to as high as 10 times the blackbody flux at 1 keV. Miller (1992) 
presented only models for T = 1 X 10 6 K, but it is likely that the deviations would be even 
larger at lower temperatures. It is not clear whether the observed X-ray spectrum and pulse 
behavior can be explained by these models, and it appears that only the hydrogen atmosphere 
is a candidate since all the heavier element atmospheres produce large flux deficits around 1 keV 
(Miller 1992). Although a two-component model may remain a necessity, it would nevertheless 
be very be useful to investigate the effects of more realistic model atmospheres on the parameters 
derived from the soft X-ray spectrum. 

5.8. Absence of Radio Emission 

Geminga has not been detected as a radio source. Models for its 180° double-pulsed 7-ray 
emission suggest that it has a much more inclined dipole than most canonical radio pulsars, 
together with a favorably inclined line of sight to the observer. With such a geometry, it is 
difficult to see why Geminga’s radio beams should be missed by observers who can see its 
7-rays. A radio-quiet Geminga may, on the other hand, be just what is expected for such a 
neutron star. Radio emission is generally thought to be powered by an inner-magnetosphere 
accelerator, perhaps even one which forms just above the polar caps where the star cannot 
supply the Goldreich- Julian charge density, p = ft • B/2nc y needed to maintain E • B = 0. 
Near Geminga’s surface, this requires a charge density n^5x 10 11 cm -3 and a possible flux 
cn ^ 1.4 x 10 22 cm -2 s" 1 . If much more charge and e^ flow than this is introduced into the 
accelerator region, they will quench that accelerator’s E-B and, presumably, the radio emission 
which it supports. 

Geminga appears to be an exceptionally good candidate for a pulsar in which 7-ray emission 
from an outer- magnetosphere accelerator produces enough e^ pairs in the inner magnetosphere 
to quench a radio-emission accelerator which would otherwise form there. This could be accom- 
plished by the e* pair production discussed in § 5.4. Intense 7-ray beams, e.g. beams 2 and 2' of 
Figure 9, also make pairs on the upward curving open field lines on which inner-magnetosphere 
accelerators are expected. A large fraction of these pairs flow in toward the star and through 
any near polar-cap accelerator. That accelerator would be quenched if the inward flowing e ± 
production rate is N± > 4 tt R 2 cn(R) ^ 2 x 10 35 s" 1 . A rate exceeding this is quite plausible 
for Geminga because of the large intensity of its inward directed 7-ray beams (e.g., beams 2 
and 2' of Figure 9). For this radio-suppressing mechanism to be effective, both a high efficiency 
for converting spindown power into 7-rays, and a large inclination so that a significant fraction 
of these 7-rays pass close enough to the star to be converted into e^ pairs, are needed. 

6. UNIDENTIFIED GALACTIC 7-RAY SOURCES 

Two recent developments in our understanding of the population of high-energy 7-ray 
sources motivate a reevaluation of the likelihood that isolated pulsars can account for the 
majority of the still unidentified COS B sources in the Galactic plane. The first is the discovery 
that Geminga and PSR 1706-44 (Thompson et al. 1992), which are older and spin less rapidly 
than than the Vela pulsar, also have higher 7-ray efficiencies. The efficiency of PSR 1706-44 
is 7% assuming isotropic emission, and it was argued here that Geminga’s efficiency is close to 
unity. The second development is the reanalysis of the COS B catalog by Mayer-Hassel wander 



& Simpson (1990) with a better diffuse background model, which leaves fewer still unidentified 
strong sources and a larger scale height for them. These trends tend to expand the pool of 
potential 7-ray sources, and to decrease the distance and hence the luminosity at which they 
are detected. 

The number of “confirmed” COS B sources above a threshold of approximately 3 x 10 “ 10 
ergs cm” 2 s” 1 is only eight. The average height above the Galactic plane for these is 2.°7, which 
is roughly twice that of the 24 sources in the 2CG catalog (Swanenburg et al. 1981). Many 
of the former “sources” closest to the plane were in reality enhancements in the diffuse flux. 
But ten new sources were suggested by the reanalysis. These have lower fluxes, and an even 
larger scale height of ~ 6°. Presumably the better background model enables fainter sources 
to be found only slightly away from the plane, so that their latitude distribution probably 
overestimates the true scale height of the Galactic 7-ray source population. If a population of 
young pulsars is to account for the COS B sources, they should have a scale height of ~ 80 pc, 
reflecting the location of their birth in a young stellar population. Adopting a scale height of 3° 
implies a typical distance of 1.5 kpc, rather than the 2-7 kpc originally assumed by Swanenburg 
et al. (1981). 

The number of potential sources with flux greater than F m j n ~ 3 X 10 10 ergs cm 2 s 1 
can be estimated from the birth rate of pulsars in the solar neighborhood, ~ 2.2 x 10” 5 yr” 1 
kpc -2 (Lyne & Smith 1990). If the apparent 7-ray luminosity is a fraction 7/ of the spindown 
power, then L^ — 4 x 10 31 t}B\ 2 P~ A ergs s” 1 , where Bu is the surface field in units of 10 12 G. 
(Any increase in apparent luminosity due to beaming is considered here to be included in 77.) 
The characteristic age of a pulsar is r = 1.6 x 10 7 P 2 B^ 2 Y r > 80 the maximum distance to which 
a pulsar of age r could be detected is r max = (L>y/4TF m \ n ) l f 2 = 5.3 77 1 / 2 (Buts)" 1 kpc, where 
75 is the age in units of 10 5 yr. The efficiency 77 is very likely a function of age, as indicated 
by the trend among the four known 7-ray pulsars. A typical value for 77 of 0.1 will first be 
examined, and then a parameterization of the form 77 = 0.2 rs will be used. 

Since pulsars younger than age 7$ would be detectable at the distance r mAX , the surface 
density of detectable pulsars at a distance r is <r(r) = 2.2 75/? = 11.6 77 1 / 2 (5{B\2 T )~ l kpc” 2 , 
where /? is the fraction of 7-ray pulsars whose beaming allows detection by us. For the case of 
constant 77, the total number of detectable pulsars to a distance D is then 

N = f 2 Tcrdra(r) = 73 ^ ^ 1 9, (11) 

Jo ^12 


where D is in kpc. A typical value of D can be the distance at which the flux of a Vela-like 
pulsar would equal F m \ ny or ^ 3 kpc. For r)(3 2 ~ 0.1 and B\ 2 ~ 3, we find N ~ 23, which can 
account for the COS B sources. 

For a model in which 77 = 0.2 75, the maximum distance to which a pulsar is detected 

can be written 
distance r 


< muuta ill wuiLll if — u.4 75, me lUOAim mu vxio laiitc muvu a puuax 10 ucvevveu 

itten as r max = 2.4 BT? r kpc, and the surface density of detectable pulsars at a 
is <7(r) = 12(Bi2r)~*f3 kpc” 2 . The total number of detectable pulsars is then 



Choosing B 12 ** 3, /? = 1, and f? m * x /Anm ~ 15 to represent the expected range of distance 
between the nearest source (Geminga) and the furthest, we find N ~ 23 as before. 

These rough calculations show that if 7-ray efficiencies of pulsars in the period range 
0.1 — 0.24 s are typical of the ones so far detected, then the standard pulsar birth rate can 
account for the 7-ray source population. This birth rate assumes a beaming factor of 5 for the 
rarfto emission. The estimates for 77 used in the calculation of the detectability of the 7- rays do 
not assume any correlation between radio and 7-ray beaming patterns, and thus do not require 
the COS B sources to be known radio pulsars. But should we worry that the detected 7-ray 
pulsars might have exceptionally high 7-ray efficiencies that are atypical of their contemporaries? 


Although this is a concern, there seems to be no strong evidence that it is the case. Pulsars of 
the age of PSR 1706-44 or younger (17,000 yr) should be distributed with a mean separation 
of 1.6 kpc. Its actual distance from us of 1.4 kpc is in line with this expectation. Similarly, 
pulsars which are younger than Geminga (340,000 yr) have a mean separation of ~ 370 pc. The 
estimated distance to Geminga of 250 pc is, therefore, also typical. Vela itself, at a distance 
of 500 pc, is closer by a factor of 4 than the mean separation of 11,000 year old pulsars. So 
the evidence is that the detected 7-ray pulsars are not overachievers, but are rather typical 
representatives of their age groups. 


7. COMMENTS 

The ROSAT X-ray observation of Geminga can be interpreted in terms of thermal emission 
from the surface of the neutron star. A combination of constraints based on fits to the X-ray 
spectrum, the brightness of the optical counterpart, and the local interstellar hydrogen column 
can be used to restrict the likely distance to the range 150-400 pc. The surface temperature is 
Ti = (5.2 ± 1.0) X 10 5 K, which is in the range of model cooling curves for a neutron star of age 
3 X 10 5 yr. The bolometric soft X-ray luminosity is 5.4 X 10 3 * (d/250 pc) 2 ergs s" 1 , but could 
be in error by a factor of 2 because of the effect of the uncertainty in Nu on the X-ray spectral 
fits. 

The implied efficiency for converting spindown power into 7-ray luminosity is then greater 
than 25%, and a substantial beaming factor could be required as well if the upper range of the 
distance estimate applies. Rather than all 7-ray pulsars having efficiencies comparable to Vela 
(tf ~ 1.4%), the observations of both PSR 1706-44 (1/ ~ 7%) and Geminga strongly suggest 
that the efficiency is a rapidly increasing function of period. The effect of this trend is to boost 
the likelihood that most of the high-energy 7-ray sources in the Galaxy are pulsars with periods 
in the range 0.1-0.24 s. They are not necessarily known radio pulsars, both because 7-ray fan 
beams can be more widely visible than narrow radio beams, and because efficient inclined 7-ray 
pulsars may quench their radio emission accelerators (§ 5.8). 

A highly inclined rotator is favored in order to account for the large 7-ray efficiency and 
the 180° separation of the 7-ray pulses. Since the 7-ray efficiency of Geminga is close to unity, 
it is of particular interest to ask how much additional luminosity might be hidden in in the 
energy band between the 1 keV X-rays and the 100 MeV 7-rays. A power law connecting the 
detections at these energies would contain only ~ 24% of the COS B luminosity, and would 
fall below the detection thresholds of all existing instruments. The claimed 2.8 a detection 
of pulsed emission from Geminga by the Figaro experiment (Massaro et al. 1993) lies above 
this interpolation, and amounts to a flux of 1.1 x 10" 10 ergs cm" 2 s" 1 in the 0.15-0.48 MeV 
range. This is only 4% of the COS B flux, but, if real, it could be indicative of an energetically 
significant component if it extrapolates over a wider energy band. It could also contain more 
than 10 38 s" 1 of e* annihilation 7-rays which may come from the pair creation mechanisms of 
§ 5.4. 

Geminga is of great importance to the study of neutron star photospheres, since it presents 
the most detailed pulse profile and spectrum so far obtained that can be confidently attributed 
to surface emission. But the X-ray pulse shapes may be difficult to analyze quantitatively. The 
single-peaked nature of the soft and hard X-ray pulses and their ~ 105° separation suggests 
a surface magnetic field geometry which includes an off-center dipole. It is not clear how 
much of the high-frequency structure in the soft X-ray pulse can be understood as rotational 
modulation of surface field, and to what extent resonant scattering in the magnetosphere may 
be responsible. The interpretation of the hard X-ray component ( T2 ~ 3 X 10 6 K) as polar 
cap emission may yet be problematic: its luminosity of 2 X 10 30 ergs s" 1 falls two orders of 
magnitude below the expected polar-cap heating from the magnetospheric accelerator. Models 
in which the polar-cap emission can be redistributed over the surface, as in § 5.4, seem to 
deserve detailed investigation. High-resolution spectroscopy and better signal-to-noise ratio in 
the harder X-rays would be a blessing (or a curse) on efforts to understand the effects of the 



surface magnetic field and any scattering processes on the emergent X-ray spectrum and pulse 
profile. 
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TABLE 1 

GEMINGA PULSE EPHEMERIS 


Parameter 

Value 

To (JD)“ 

2448329.886823149 

P 00 

0.2370973861 

P (s s *) 

1.0977 x 10~ 14 


4 TDB time of phase 0 in Figures 1 and 2. 


TABLE 2 

FIT TO DOUBLE BLACKBODY MODEL 


Parameter 

Value 

Confidence 

T\ 

5.2 (+1.0-1.3) x 10 s K 

68% 


(+1.2, -2.2) 

90% 

t 2 

3.0 (+0.7, -0.6) x 10 6 K 

68% 


(+0.9, -0.7) 

90% 

N H 

1.5 (+0.8, -0.5) x lO 20 cm" 2 

68% 


(+1.7, -0.6) 

90% 

/ 

3.5 (+4.5, -2.9) x lO" 5 

68% 


(+6.6, -3.4) 

90% 

C 

162. 


Fx 

7.2 x 10 -12 ergs cm -2 s -1 


Ft 

2.6 X 10 -13 ergs cm -2 s -1 



4 Bolometric flux in component before interstellar absorption. 


TABLE 3 

FIT TO BLACKBODY PLUS POWER-LAW MODEL 


Parameter 

Value 

Confidence 

T i 

4.5 (+1.4, -1.5) x 10 5 K 

68 % 


(+ 1 . 6 , - 2 . 1 ) 

90% 

a 

1.47 (+0.28,-0.37) 

68 % 


(+0.38,-0.55) 

90% 

N h 

1.85 (+1.3, -0.7) xlO 20 cm " 2 

68 % 


(+2.4, -0.8) 

90% 

f 

1.93 (+9.1-1.9) xlO -7 

68 % 


(+12.1-1.9) 

90% 

C 

545. 


F? 

1.1 x 10 -11 ergs cm -2 s _1 



* Bolometric flux in blackbody component, before interstellar absorption. 


TABLE 4 

SUMMARY OF ENERGY BUDGET 


Component 

Luminosity (ergs s 1 ) 

Spindown Power 

3.3 x 10 34 

7 -rays 0 

2.3 x 10 34 

“Soft” X-rays'* 

5.4 x 10 31 

“Hard” X-rays'* 

1.9 x 10 30 


* Luminosity for a distance of 250 pc, in the isotropic case. 
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FIGURE CAPTIONS 


FIG. 1.- X-ray pulse profile folded according to the ephemeris of Table 1. The top panel shows 
all the data in the 0.07-1.5 keV range, and the lower three panels show the same data 
separated into three energy bands. The data are repeated for two cycles to guide the eye. 
Phase 0.0 corresponds to To from Table 1. 

FIG. 2.- A comparison of the X-ray and 7-ray pulses in absolute phase. Phase 0.0 corresponds 
to To from Table 1. The 7-ray light curve is taken from Bertsch et al. (1992). 

FIG. 3.- Fit of the double blackbody model to the X-ray spectrum of Geminga. (a) Raw counts 
and model spectra folded through the detector response matrix, (b) Unfolded spectral 
model. 

FIG. 4.- Confidence contours for the parameters of the soft X-ray component in the double 
blackbody fit. The confidence levels correspond to the case of three interesting parameters 
(Ti , A'h, and C). Visual magnitudes are extrapolated from the Rayleigh- Jeans tail of the 
blackbody fit, while distances assume emission from the full surface of a neutron star of 
radius R = 10 km. 

FIG. 5.- Confidence contours for the parameters of the hard X-ray component in the double 
blackbody fit. The confidence levels correspond to the case of two interesting parameters 
(T2 and /). Corresponding values of T\ are given as dashed lines. 

FIG. 6.- Fit of the blackbody plus power-law model to the X-ray spectrum of Geminga. (a) 
Raw counts and model spectra folded through the detector response matrix, (b) Unfolded 
spectral model. 

FIG. 7.- Confidence contours for the parameters of the soft X-ray component in the blackbody 
plus power-law fit. The confidence levels correspond to the case of three interesting 
parameters (Ti,iV H , and C). Visual magnitudes are extrapolated from the Rayleigh- 
Jeans tail of the blackbody fit, while distances assume emission from the full surface of a 
neutron star of radius R = 10 km. 

FIG. 8.- Proposed geometry of the outer magnetosphere of Geminga, showing the accelerators 
and 7-ray beams powered by them. The accelerator, which is relatively empty of charge, is 
the stippled region. The accelerator on the right (left) is bounded from below (above) by 
the last closed field-line surface. The open field-line bundle between the accelerator and 
the polar cap (vertically shaded region) contains N, the e _ (e + ) flow out of the accelerator 
down toward the polar cap (see also Figure 9). 

FIG. 9.- (a) Geometry of the surface magnetic field, showing the highly inclined sunspot dipole, 
and 7-ray and particle beams inwardly directed along the open field-line bundle, (b) Top 
view of the neutron star, showing the near-surface field configuration and the rotation 
phase for which an observer at the bottom of the figure would see the maximum hard 
X-ray flux. The dashed lines are the radial directions in which | E x B | can vanish at 
r ~ 3 R. 

FIG. 10.- Schematic phase diagram of the 7-ray and X-ray light curves according to the theory 
presented in §§ 5. 1-5.5. Qualitative similarities with the observational data shown in 
Figure 2 are apparent. 
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